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ABSTRACT 


Few underwater acoustic models have been validated in the Very Low Frequency (VLF) 
regime. The purpose here is to explore one model’s abilities and compare it with ocean data once 
it is ascertained the model has converged to its final solution. A comparative analysis of data 
recorded in NATIVE 1 (the Noise and Transmission Loss in VLF Environments) Event 1 and results 

: generated using the Seismo-Acoustic Fast field Algorithm for Range-Independent Environments 


(SAFARI) is the focus of this study, with emphasis on understanding some of the model’s features. 
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I. INTRODUCTION 

In submarine warfare there is a need for accurate acoustic 
models both for the prediction of how well a given sensor will 
work and to support the design of new sensors. These acoustic 
models need to be verified, in the acoustic regime of 
interest, against a database of experimentally-measured data 
representative of the environment. | 

The verification of an acoustic model requires three 
primary ingredients. The first is a model that produces 
results that are testable, the second is a database against 
which to compare the results of the model, and third is a 
method of comparison; whether statistical or graphical. One 
purpose of any acoustic model is to simulate the environment 
with sufficient accuracy so results can be compared to 
measurements. Limitations on the model input parameters have 
direct ramifications on its results. Similarly the precision 
and accuracy of the equipment used in recording the data as 
well as any perturbations which may occur during the 
collection of the data must be taken into account when 
determining its quality. A model is only as good as its 
theoretical design and an experiment is only as good as re 
design and equipment. Once both results are compiled they can 


be compared. Then a determination can be made as to whether 


the model matches with the features of the measured data. 





The objective of this research is to verify the Seismo 
Acoustic Fast field Algorithm for Range-Independent 
environments (SAFARI) in the very-low-frequency (VLF) acoustic 
regime (from 5 Hz to 20 Hz) by comparison with the data 
recorded by Marine Physical Laboratory (MPL) in the Noise and 
Transmission Loss in Very Low Frequency Environment (NATIVE 1) 
experiment, which they conducted in August of 1990. Chapter 
II gives a thorough background on the NATIVE 1 experiment, 
including a description of the site’s geoacoustic and 
bathymetric makeup as well as the types of source and receiver 
used. Chapter III goes into a description of the theory 
behind SAFARI. Chapter IV discusses the data sources used in 
the comparison. Chapter V investigates model aliasing and 
contour offset features. Chapter VI presents a graphical 
comparison of model results and data. Chapter VII gives the 


conclusions reached and proposes areas for further 


investigation. 





II. NATIVE 1 


A. BACKGROUND 

In August of 1990, the NATIVE 1 experiments were conducted 
to acquire data to provide insight into VLF-propagation 
effects. Along with the development of a VLF data base was 
the effort to design VLF performance prediction models which 
would support the development of new Antisubmarine Warfare 
(ASW) acoustic sensors. NATIVE 1 was a multi-laboratory,at- 
sea exercise conducted on Blakes Plateau, Figure 1 
(Ref. 1}. Two sites were used, one in deep water and 
one in shallow water. Five events, lasting several weeks were 
conducted using a variety of devices and equipment. Data from 
the first of these events at the deep water site will be used. 
Figure 2 (Ref. 2] shows the locations of the initial 
entry and recovery points of the receivers. Typical Swallow 
Float deployments consist of forming an equilateral triangle 
_ pattern with three floats, at known positions, spaced about 
6.5 km between floats. This provides a reference for float 
localization. A 15-lb anchor on a 3-m tether replaces the 
ballast. The remaining floats are deployed in the center of 


the triangle. Figure 3 shows the path of the source over the 


day and a half comprising the recording time. The event 
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Figure 1 Native 1 Site "D" in 3400 m of water and "S" in 
1000 m of water.[(Ref. l:p. 2] , 


starts with the source well north, approximately 12 km, of the 
center of the sensor site. It maintains a track of 
approximately 193 degrees past the receivers and continues 
south for another 50 km, and then turns around on the 
reciprocal heading finishing approximately 6 km north of the 
site. 
B. EXTREMELY LOW FREQUENCY SOURCE 

In this event a calibrated extremely low frequency (ELF) 


source projector was towed by the research vessel (R/V) Gyre. 


A precise positional log (Appendix A) was maintained from 
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Figure 2 Plan view of entry and recovery positions of 
Swallow floats. Event 1, NATIVE 1, (34°N, 74°W) [Ref 2:p. 2] 
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Figure 3 Research vessel Gyre track during Event 1, NATIVE 
1 (31°N, 74°W) [Ref. 2:p. 2] 

multiple navigational sources. These included Global 
Positioning Systems (GPS), LORAN, and Transit Satellite 
Navigation (Transit Sat Nav). Towing speed was 4.0 + 0.5 


knots and the source depth was 122 + 7 meters. The source 


projected two continuous wave frequencies simultaneously; 7 Hz 





and 10 Hz from 1630 to 2230 on 4 August and 0030 to 0400 on 5 
August, and 10 Hz and 16 Hz from 0800 to 1130 on 5 August. At 
7 Hz, the average source level was 164 + 2 dB//uPa, at 10 Hz 
the average source level was 166 + 2 GB/uPa, and at 16 Hz the 
average source level was 161 + 2 adB/uPa (Ref. 3]. 
(All estimates were made from visual inspection of the data.) 
The positional accuracy of the ship was source dependent with 
the GPS system having the highest degree of accuracy, about + 
60 meters. Transit Sat Nav was next with an accuracy of + 100 
meters and LORAN was accurate to within + 0.5 nautical miles. 
The times each of these systems were used on the above dates 
are highlighted in Appendix A, along with the actual track 
data. 


"The Transit Sat Nav is dependent on satellite 
elevation above the horizon, time from satellite fix 
to position solution and dead reckoning parameters 
(course and speed) used for updates. Therefore I 
estimate Transit Sat Nav is only accurate in general, 
to within 0.5-1.0 nautical miles. LORAN appears to be 
accurate to within 0.5 nautical miles with better 
repeatability. For this reason there is a constant 
offset or position error between GPS and LORAN of - 
0.021 degrees longitude which should be added to LORAN 
longitudes. This would move all LORAN positions 
further to the west. This correction has not been 
applied to LORAN positions in enclosure(2) (Appendix 
A) which is taken from the NADC navigation 
rectification products. It also will not affect 
ranges along the track from source to receiver due to 
the north-south orientation of the Event 1 track." 
(Ref. 4] 








C. ENVIRONMENTAL DATA 

The exercise area is located off Jacksonville Florida, 
east of the Gulf Stream in the waters commonly called the 
Sargasso Sea. At the end of Event 1, 41 bathythermographs and 
one conductivity, depth, temperature (CDT) cast were made 
(Figure 4). Sample values of those results [Ref. l:pp. Bi- 


B23] are in Appendix B, and are typical of the water masses 
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Figure 4 Site "D" temperature profiles for (a)AXBTs, (b) XBTs 
[Ref. 1l:p. 22] 


variability was minimal with one exception located in 
Gulfstream Waters well to the west of the track. Province 
contours (Figure 5) were generated and representative values 


of each province are listed in Appendix B. Those profiles 


PATIVE 7 SOUHO SPEED PROFILE PROVINCES 





LONGITUDE (DFG) 
Figure 5 Site "D" sound speed provinces.[(Ref. l:p. 23] a 


that are generally cooler have lower province numbers. Table 
I shows the province profile number and range from site "D" 
along the track of R/V Gyre during Event 1. 

Seasonal heating in summertime dominates the surface water 
while the seasonal thermocline extends to about 
100 m. From 100 to 500 m is the regime called ‘18 degree’ 
water mass where the temperature gradient is small producing 
an isovelocity water. Below 500 m to any where between 2500 


to 3000 m North Atlantic Central Water with a deep sound 


Table I RANGE FROM SITE “D"" AT WHICH THE SOURCE ENTERED A 
PROVINCE PROFILE NUMBER ALONG EVENT 1 TRACK. [Ref 1:p.24] 


Range (nautical miles) 
12.54 ae ees 
14.21 















channel axis at about 1200 meters. Below this is the very 
stable North Atlantic Deep and Bottom water mass. 
D. GEOACOUSTIC BOTTOM 


The bathymetric track (Figure 6) was recorded by R/V Gyre 
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Figure 6 Bathymetric Track from Site "D", Event 1 [Ref. 1:p. 
40] 


during Event 1. The bottom varied from a minimum of 3159 m 
near the VEDABS 1 location (center of the site) to a maximun 
at 3688 m. This yields an approximate slope from 0.002 to 
0.001 for the bottom and, as such, is considered horizontal 
and flat. In Figure 7 representative samples of the 
geoacoustics of the bottom sediments from ([(Ref. 5] 
were generated along the R/V Gyre track using the model found 
in the front of Appendix C, [Ref 1:p. 45] 
E. SWALLOW FLOATS 

The following discussion is adapted from D’Spain et.al. 
[Ref. 6]. A general description of the Swallow Float 
is followed by a discourse on the data acquisition process, 


including a discussion of the hydrophone and the electronic 
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Figure 7 Site "D" Bearing OOO°T Water Depth and Depth to 
Acoustic Basement[Ref. 5:p. B-24] 


circuit, and finally a discussion on the hydrophones 
autospectra. 
1. Description 

Swallow Floats are independent, freely-drifting, 
neutrally-buoyant sensors which can be ballasted to float at 
any desired depth (+ 100 m). Figure 8 shows a schematic of a 
typical Swallow Float. Mounted within the 0.432m diameter 
Benthos glass shell are three orthogonally mounted geophones, 
a compass, the electronics with power supply (not shown), and 
cassette tape recorder (not shown). A polyethylene hard hat 
covers the glass shell’s exterior with an Ocean Applied 
Research (OAR) xenon flasher attached to one side and an OAR 
radio beacon on the other. These are used for detection and 


localization during recovery. 
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Figure 8 Swallow float TRef. 6: ip. 31] 

To prevent distortion of the compass readings, the 
compass is mounted on the bottom of the shell at its maximun 
distance from the magnetic fields of the geophones, mounted on 
the top of the shell. The power supply provides a maximum of 
200 hours of data acquisition and 25 days of aeoustie: 
transponder life. Duration of Swallow float deployments is 


not limited by power supply requirements but by the capacity 


of the 17-megabyte cassette tape recorder. 





An infrasonic hydrophone is mounted outside the glass 
shell on the base plate just below the compass, approximately 
0.5 m below the geophones. Suspended below the base plate is 
the expendable ballast. The ballast can be acoustically 
commanded to release and has two independent backup timed 
release circuits. The ballast weight is accurately determined 
prior to float deployment to achieve the desired depth of 
deployment. A change in ballast of 1 gm results in depth 
changes of 7 to 10 m. The ballast fits snugly against the 
base plate in order to reduce its rocking motion relative to 
the float; herein Figure 8 is distorted. 

Each float is also equipped with an International 
Transducer Corporation (ITC) localization hydrophone which is 
suspended 1.83 m below the glass shell. This hydrophone has 
a source level of 191 dB//luPa at 1 m and a receiving 
sensitivity of ~180 dB//1V/uPa at 8 kHz. Under microprocessor 
control and synchronized prior to deployment, the floats send 
8-kHz localization pulses of 10-ms duration. This allows for 
post deployment determination of depth and location during 
deployment. This hydrophone and the associated circuitry also 
allow for two types of acoustic commands, uniquely coded for 
each float, to be issued. The first is a 50-ms 8-kHz pulse 
used for ship to float slant range calculation and the second 
activates the ballast release circuitry and issues a 500-ms 8- 


kHz pulse and is used at the end of the experiment for float 


recall. 








2. Data Acquisition 
The Swallow Float data acquisition system (Figure 9) 
records the acoustic pressure and particle motion at the 


sensors. This discussion is restricted to the acoustic 


BLOCK DIAGRAM OF THE SWALLOW FLOAT 
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Figure 9 Swallow Float data acquisition system [Ref. 6:p. 
31] 
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pressure only. An OAS model E-4SD hydrophone is the pressure 
sensor which converts acoustic and nonacoustic pressure 
fluctuations into voltage signals with a sensitivity of -182 
aB// 1V/upPa. The hydrophone has a frequency response which 
is flat (+ 1 dB) from 0 to 5kHz [Ref 6:p. 8]. The signal then 


passes in order through: 


e two prewhitening RC filters 
e fixed 80 dB gain amplifier 


@ three highpass RC filters and four lowpass RC filters 
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e automatic gain control (AGC) amplifier 
e anti-aliasing filter 

e A/D converter | 

edata buffer with 44 seconds of data 


® cassette tape recorder with 1 second of write time for 
each acoustic record 


Each acoustic record is 45 s long (44 seconds of recorded 
data and 1 second of write time). 
3. Pressure Autospectra 
From the recorded data, the RMS pressure at each float 
was calculated. These RMS time series were corrected for AGC 
gain but were otherwise uncalibrated. The vertical axis of 
the recorded data is in units of RMS volts at the A/D 
converter (Figure 9). 


"Each RMS level value results from taking the square root 
of the average of the squared amplitude levels over a 
period of five seconds (or 250 points since the data 
sampling rate is 50 Hz). Since each record is 2250 points 
(45 seconds) long, then nine RMS power values are 
calculated for each record. The RMS value for the last 
five-second period of each record includes one second of 
zeros. The one second of zeros represents the time during 
which no data is sampled while data in a temporary buffer 
is being written to cassette tape." [Ref. 2:p. 21] 


"Spectral estimates with two different frequency 
resolutions are calculated. Those with a 98 mHz 
resolution are made by dividing 40.96 seconds in a data 
record (obtained after skipping the first 3 seconds of 
data) into seven 10.24-second-long segments, with a 50 
percent overlap between segments. The segments are fast 
Fourier transformed after using a Kaiser-Bessel window of 
a=2.5 to window the data. The resulting seven spectra 
within a record are incoherently averaged together. The 
averaged spectra for four consecutive records (equal to 
three minutes of data) are then incoherently averaged in 


14 


order to further reduce the variance of the spectral 
estimates. For spectral estimates with a 24 mHz 
resolution , the segment length which is transformed is 
quadrupled to 40.96 seconds and the spectra for twenty 
consecutive records are incoherently averaged."[{Ref. 2:p. 
26) 


The autospectra are calculated in dB // 1(“Pa)2/Hz. Records 
used in the analysis were all three minute averaged data, 


using four 45 second records for each data point. 
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IIIT. MATHEMATICAL MODEL 


A. BACKGROUND 

In the attempt to model the complex environment found in 
the ocean, there have been a variety of techniques used to 
determine the solution to the wave equation, each involving 
its own simplifying assumptions and approximations. One 
subset of these techniques is the method using separation of 
variables. Two approaches fall into this general category: 
the normal mode and the fast field approach. It is the latter 
that is implemented in SAFARI ([Ref. 7). The following 
sections provide a summary of the wavefield technique found in 
Schmidt [Ref 7}. In a horizontally stratified environment, 
including fluid and solid layers, this technique yields, an 


exact solution to the linearized wave equation, 
ees eel 
e(r,z)? dt? 


)P(zr,z,t)=F,(r,Z, t) ee 


The following sections develop the basic mathematical 
approach, discuss environmental and analytical models, 
boundary conditions, and numerical techniques including, the 
global matrix approach, wavenumber integration, aliasing, and 
contour offsetting. 
B. MATHEMATICAL MODEL 

As early as 1948 [Ref. 8}, the subject of wave 


propagation in a stratified medium was being studied in planar 
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two- and three-layer environmental models. The technique 
applies a series of integral transforms to the linearized 
Helmholtz wave equation, Eq. 3.1, reducing a partial 
differential equation in time and three spatial dimensions 
into a series of ordinary differential equations in the depth 
coordinate. These ordinary differential equations can be 
solved analytically in terms of unknown amplitudes by matching 
boundary conditions at the interfaces. For a cylindrical 
coordinate system {r,8@,z}, the model solves for an axial 
symmetric, or planar, wavefield independent of ee. The 
isotropic media can be expressed in terms of scalar potentials 
V(r,z,t) which satisfy Equation 3.1, where c(r,z) is the wave 
speed and F,(r,z,t) is a forcing term representing the source. 
The first of two integral transforms is next applied. The 


forward Fourier transform 


tf -iewt 
f(o) ig [Fieve dt (a) 


3.2 
F(t) =f £(o) e#**at (b) 
gives the frequency domain wave equation 
(W+ka(r,z)) ¥(r,z,) =f,(r,z,0) 3.3 


where w is the angular frequency and k,(r,z) is the medium 


wave number, 


k, (2,2) =0/c(r,z) 3.4 
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For stratified range independent environments with a point 


source present, Equation 3.3 can be simplified to 


£,(z,@) 6 (r) 
2nur 


(W+ki(z)) ¥(r, 2,0) = 3.5 


In the following, the w dependence will be suppressed. The 


Hankel transform pair are utilized 


g(k) =[G(r) J,(kr) rdr (a) 
0 


G(r) =fo(x) J,(kr) kdk (b) 
0 


where k is the horizontal wavenumber, and the forward 
transform, Eq 3.6(a), is applied to obtain the depth-separated 
wave equation which is an ordinary differential equation in 
depth only. The solution to this equation is the sum of a 
particular solution ¥, to Eq. 3.5, and any linear combination 
of the two independent solutions to the homogeneous equation. 
It is commonly called the depth-dependent Green’s function 
whose coefficients are arbitrary and are determined by the 
boundary conditions. Once the coefficients are found, the 
inverse Hankel transform (Eq. 3.6(b)] yields a wavefield at 
any range r at the angular frequency w. The inverse Fourier 


transform yields the time response in the same fashion. 
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C. BOUNDARY CONDITIONS 


1. Environmental Model 


In SAFARI, the environment 


is represented as a 


horizontally stratified medium, Figure 10, with all interfaces 


planar and parallel. 
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Figure 10 Horizontally stratified environment [Ref. 7:p. 5] 


independent. 
isotropic with Lamé constants \ and 


compressional wave speed is 


co=¥ (A+2p) Tp 7 


and the shear speed is 


cy=VE/p 


Solid layers are assumed to be homogeneous and 


&@ and density p. The 


Physically realistic materials must satisfy the condition 


c,<¥0.75c, 


3.9 





so that the bulk modulus, K= A+2/3u, be positive, and the 
solid medium have positive compressibility. The density in 
fluid layers must be constant while the nonvanishing Lamé 


constant \ is allowed to vary with depth according to: 


1. 
atone 3.10 


corresponding to the following depth dependence for the sound 


speed, Eq 3.7, 


C,(2) =/17 (p (az+5) ) 3.11 


2. Analytical Models 

The analytical field representation for the media 
included in SAFARI are for a homogeneous fluid, inhomogeneous 
fluid, and solid representations, and are governed by the 
equation of motion for each media. The focus of this section 
is not on the details of each derivation but rather on the 
fact that the wavefield at each interface has two integral 
representations one each from the upper and lower layers. 
For the homogenous fluid, the solution to the depth- 


dependent Greens function for the potential is of the form 


@ (r, 2) =[ [A"e**+A*e*tU, (kr) kdk 3.12 
0 
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where k,(z)=h, (a constant), &(k)=/k?-hn , and no source is 


present. A’ and At are the unknown amplitudes of the downgoing 
and upgoing waves. If an omnidirectional point source is 
present the particular solution is added to the homogeneous 


solution and is of the form 


Ss BA e @l2-Z6l 
a Oe) 3.13 
® (r,z) oa —— Jo (kr) kdk 


where S, is the source strength and z, is the source depth. The 


associated displacement and stress are 


S27 Jie 
w(r,Z) =- 2 [sign(z-z,) e “l?-2l 7, (kr) kdk 3.14 
ie] 


-a|z-z 
2 ee — Jq (kr) kdk 3.15 


0,,(r,2) =-—-——_ 


where 


sisate-2) + 


Since the ocean is not a homogeneous media and the 
sound velocity varies with depth, a single. layer ocean wave 
guide doesn’t properly reflect acoustic propagation. However, 


if a range independent wave guide is represented by an 
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number of homogeneous layers, a numerical solution based on 
their field representation will converge toward the correct 
solution. It can be shown when modelling a deep water site 
there must be a velocity profile incorporated, as is the case 
here, moreover, to ensure a satisfactory convergence the 
isovelocity layers must be less than one quarter wavelength 
thick. Once again it can be shown the solution to the depth 
dependent Greens function, based on c(z), from Eq. 3.11, and 
not an isovelocity constant, can be obtained by using a 
variable transformation in which 
C=c 2/3 (cz+d) 3.16 


= (pw?a) ~2/3 (k?2-pw? (az+b) ) 


When substituted into Eq. 3.5 the following is obtained: 


(2-7) (7) =0 3.17 
“dt? ; 


This yields the integral form of the potential with a source, 


aie Na So y 
® (r,z) = Gao (ke) kek 


2073 (Ai (€,) +iBi(@,) ) Ai (2) 
Ai!(€,) (Ai(€,) +iBi((,)) -Ai(¢,) (Ai/(€,) +iBi’(¢,) ) 


,a(Zz-zZ,) <0, 


3.18 


2e7°V3Ai (C,) (Ai (f) +iBi(f)) 
Ai‘((,) (Ai (€,) +iBi((,)) -Ai(€,) (Ai/(@,) +iBi’(¢,) ) 


,a(Z-Z,) 20 
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By expressing the displacements due to both normal 
stress and tangential stress forces in terms of scalar 
potentials the solution to the depth dependent Greens function 


for a solid are of the, form, 
®(r, 2) =[ [A-e-*#+A*e%*] Jy (kr) kdk 3.19 
GQ 2 


¥(r,2Z) =[ [B-eB+B*eP*) J, (kr) dk 3.20 
0 


The normal stress and tangential stress are solved for 


in the same fashion as in the homogeneous case. Here 6 is the 





wavenumber for tangential stress analogous to the normal 
stress wavenumber a in Eq. 3.12. 
3. Matching Boundaries 
The boundary conditions are matched at each interface 
dependent upon the type of interface involved. The following 
list summarizes the conditions: 

e For an interface separating two fluid layers, the vertical 
displacement w and the normal stress o, have to be 
continuous. If one of the media is a vacuum, the normal 
stress must vanish. 

os e For an interface separating a solid and a vacuum, both oa, 
and o, must vanish, whereas a solid-fluid interface 
requires that w and o, must be continuous and oa, must 
vanish. 


® For an interface separating two solid media, w,u,o,,, and 
o,, must all be continuous. 


"Since the boundary conditions have to be satisfied at all 
ranges r, it is obvious that they must be satisfied by the 
kernels in the integral representations as well. The 
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radiation conditions for z»to together with the conditions 
to be satisfied at all interfaces simultaneously, lead to 
a linear system of equations in the unknown kernel 
coefficients A,A*,B, and Bt. In principle, this system 
has to be solved for all values of the horizontal 
wavenumber k, and the total field can then be determined 
by evaluating the inverse transforms. Except for a few 
trivial cases, however, both the solution of the linear 
system of equations and the evaluation of the inverse 
transforms have to be done numerically, requiring 
truncation and discretization of the horizontal wavenumber 
axis." (Ref. 7:pp. 17-18] 


D. NUMERICAL TECHNIQUES 
The global matrix approach and wavenumber integration are 
. now discussed. 
1. Global Matrix Approach 
The global matrix approach is a more direct and 
computationally more efficient technique in assembling the 
local solution to the depth dependent Greens function, for 
each layer, in terms of its source contribution and unknown 
scalar potentials, with all other local solutions. The 
boundary conditions determine a system of equations in the 
Hankel transforms of the potentials to be satisfied at each 
interface. In similar fashion to finite element model 
programs [{Ref. 9] these local system of equations are 
mapped into a global set of equations defining the boundary 
conditions at all interfaces simultaneously. 
"The resulting global coefficient matrix is organized to 
be block-bidiagonal and diagonally dominant in close 
analogy to the global stiffness matrix arising from the 
finite- element method. Unconditionally numerically- 
stable solutions are determined efficiently by gaussian 


elimination, yielding the field in all layers 
simultaneously." [Ref 7:p. 29] 
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The mechanics of the technique’s derivation of the solution is 
secondary to the understanding that unconditionally stable 
solutions are obtained, and that the technique is stable for 
any machine precision. Herein, SAFARI is an excellent choice 
for depth range contouring of single frequency transmission 
loss, the focus of this study. 
2. Wavenumber Integration 
Numerical techniques which evaluate the depth 


dependent Green’s function by the inverse Hankel transform, 


G(r, Z) =[o(k, Z) J, (kr) kdk 3.6(b) 
Q 


derive truncated solutions to the infinite integration 
interval. Here the subscript m equals 0 except for the shear 
stress and horizontal displacement where m equals 1. Because 
the interval on the wavenumber k goes from 0 to » the solution 
can be done to any degree of accuracy within the computing 
limits of the machine. In a lossless media, waveguide 
propagation produces normal modes and interface waves which 
show up in the wavenumber integration as poles. As the 
integrand varies over k these plot in the complex plane on the 
real axis. For lossy media, these poles move out into the 
complex plane, remaining close to the real axis. In cases 
where solid layers underlie fluid layers care must be taken to 


ensure the maximum medium wavenumber truncation point is 


selected beyond any significant contributions from waves 
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propagating in those solid layers. Here, evanescent surface, 
interface or plate modes may exist causing poles to appear in 
the integration plots. Typically 20% above a maximum medium 
wavenumber is a good value for truncation of the wavenumber 
range. 

Once the truncation point is selected two factors 
complicate the remaining solution. The first is the existence Ps 
of normal modes and interface waves which generate poles in 
the kernel of the transform. The second is the Bessel 
function J,,(kr) which may vary rapidly in k, the horizontal 
wavenumber, for long ranges in r. A very fine sampling 
interval has to be chosen to avoid aliasing, addressed below. 

DiNapoli and Deavenport (Ref. 10] show, for 


the ranges of this study, the Bessel function, 


J,, (Kr) =+ (HEY (kr) +H? (kr) ) 3.21 


can be represented by the outgoing Hankel function H,,” 


-i(kr-(m++) tn) 
lim Hf? (kr) | Be 22 3.22 
ukr 


kr- 


(The incoming Hankel function is needed only for evaluation in 


a 


the near field within a few wavelengths of the source.) This 


yields an inverse Hankel transform integral 


| 1 img) ore -ik 
G(r, z) =| —— k, e ~ikrdk, 3.23 
(zr, 2) sare fo z) Vk 
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The wavenumber space having been truncated is divided 
equidistantly in SAFARI as follows: 


k,=Kyptlak, 120,1...(M-1) 3.24 


where M is the total number of sampling points. The range 
axis is divided accordingly as 


+IAE; j=0,1...(M-1) 3.25 


Ty5*Inin 


where the range step Ar is governed by the relationship 


ArAk=2n/M 3.26 


and M is an integral power of 2. The discrete approximation 


of Eq 3.23 is finally obtained, 


G(r,,z) OOP cae 9 [g(k,, 2) @HaolBk 1] 9A (2nls/im 
3.27 
The summation can be performed by means of a fast Fourier 
transform (FFT) generating the field at all M ranges, in Eq. 
3.25, simultaneously. 
3. Aliasing and Contour Offset 
Regarding aliasing, the undersampling in one domain 
causing wrap around in the other domain of a transform, Eq. 
3.27 does not yield G(r,z) but a summation from the lower 
limit to the upper limit of the interval. The range from zero 
to the lower limit and from the upper limit to infinity will 


cause some aliasing in the wavenumber. Over the interval 


[Lminr Umn + R], Where the range window width R=MAr, the maximum 
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useable range is always R=MAr, independent of the choice of 
Lin ° 
For ranges larger than r,,+R aliasing can be reduced 
by choosing the window R so large that the signal is known to 
Gissipate within the window. For a fixed wave number interval 
this would require a large number of wavenumber samples. 
According to Cauchy’s theorem the integrals in the 
complex plane between two points are invariant to a change in 
the integration contour, that is by moving the integration 
contour away from the pole, offsetting it, the integrands can 
be smoothed and aliasing avoided. In SAFARI this is done by 
replacing k with k=kt+ie where € is a vertical offset whose 


contribution to the total integral is insignificant when 


compared to that of k ie. e«k,,.-k,,. Eq. 3.27 is replaced by 


1 fim 5) ge -ik 
/ -\ i k, e-iktdk 3.28 
G(r, Z) oar’ fo Zz) VK 


where C is the contour shown in Figure 11. For ranges larger 
than r,,, + R, a factor of e“® is introduced to the FFT and the 
signals will be attenuated. Similarly for ranges smaller than 
Ynin Signals wrapped around will be amplified by e*®. No change 
occurs in the maximum range r,,,=R, from the real axis 
integration. Corresponding to an attenuation of 60 dB of the 


wrap around signal the value of e is formulated as 


E€= 





3 3 
= -k. 3.29 
Rloge 2x(M-1)loge (Kau Kin) 
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Poles 


Ree cut for Vk?-kz 


Figure 11 Complex Integration Contour [Ref. 7:p. 33] 


which even for a small number of sampling points ensures ¢€ « 
Kuax7Kmine 4Paramount to this numerical integration scheme is 
that the kernel must vary smoothly between sampling points. 
If the integration is truncated at a point where the kernels 
amplitude is significant, the discontinuity will give rise to 
aliasing. To avoid this, a first order Hermite polynomial is 
employed and the kernel and its first derivative become 


continuous at the truncation point. 
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IV. DATA SOURCES AND ACOUSTIC MODELING INPUTS 
The recorded transmission loss data sent by Marine 
Physical Laboratory (MPL) of Swallow Float 0 is described 


along with the input files used in running the SAFARI model. 


A. SWALLOW FLOAT DATA 
The data [Ref. 11] received from MPL were 12 
files of Swallow Float 0 data recorded during Event 1 of the 
NATIVE 1 experiment. It was received via electronic mail and 
transferred to 5.25" high density floppy disks. The files 
contained two columns of information. The first column 
contained the R/V Gyre-to-buoy range in kilometers and the 
second column contained either the percentage of clipped data 
points (four files) or the transmission loss (eight files) in 
units of dB referenced to a source with pressure amplitude of 
1 wPa at 1 meter. Reference 11 received from MPL goes into 
the scheme for defining the data: TLO.[1,2,3,4].[7,10,16] 
references the Transmission Loss Float 0, the region number 
and the frequency for the region. The regions are, 
Region 1: The Gyre approached the float from the north- 
edna near the beginning of Event 1 broadcasting 7 and 


Region 2: The Gyre passed over the float heading south- 
southwest broadcasting 7 and 10 Hz. 


Region 3: After the SUS charge events as the Gyre continued 
to the south-southwest broadcasting 7 and 10 Hz. 
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Region 4: The Gyre approached the float from the south- 
southwest broadcasting 10 and 16 Hz. 


With regard to the navigational data and the Gyre-to-float 
range, the 5 minute log in Appendix A was interpolated every 
45 seconds to get a position for each record of data. 

As of the writing of this paper, there is a question as to 
whether or not a range offset of approximately 0.5 km occurs 
in the navigational data. 

Uniess otherwise annotated, all Swallow Float data 
presented in the analysis section has less than 5 percent of 
clipped time series points in the hydrophone channel. Where 
the percentage of clipped points exceeded 5 percent the record 


was deleted. This was done using a DOS text editor. 


B. SAFARI MODEL AND INPUTS 

The SAFARI model was installed at the Fleet Numerical 
Oceanography Center in Monterey, California in a UNIX 
environment on an HP workstation. It was supported by DISSPLA 
graphics, and executed the Fortran code on an F77 Fortran 
compiler. Upon satisfactory duplication of the results for the 
test files in Ref. 7, and an independent two-layer split fluid 
case verified by alternate means, analysis proceeded with the 
model. | 

The water column data from province 4 in Appendix B and 
the geoacoustic model data from Appendix C was used to produce 


a 43-layer water-sediment input file in the standard SAFARI 
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FIP format. The FIP module generates single-frequency 
calculation of the total wavefield for a single point source 
whose source strength S, is normalized to yield a pressure of 
1 Pa at 1 meter distance from the source. The graphical 
results correspond to the standard definition of transmission 
loss. [Ref 7:p. 67] The table in Figure 12 shows the required 
input file blocks. The first block contains the title of the 
run, the second block the options selected. For this study 
the following options were used at one point: 


@® N the calculation of normal stress o, (equal to negative 
pressure in a fluid). 


® I Hankel transform integrands are plotted versus 
wavenumber. 


® T Transmission loss of the pressure versus range. 

e J Complex integration contour. 

@ Z Plot of the water column velocity profile. 
The third block contains the frequency which was either 7 Hz 
or 10 Hz, and the integration contour offset. If left blank no 
contour offset was incorporated. A zero input gave the 
default contour offset which yields 60 dB at the longest range 
considered in the FFT. The fourth block is the number of 
layers (NL) including the upper and lower halfspaces, and is 
followed by each layers upper boundary depth from 1 to NL, the 
compressional velocity, the shear velocity, the compressional 


attenuation in dB/A, the shear attenuation in dB/A, (both of 


which are converted from dB/m*kHz to @B/A by multiplying each 

















Block Parameter Units Limits 
errr pene ee A SS SS 
I TITLE : title of run - < 80 char. 
nn 
I opti opt2 ...: output options - < 40 char. 
er PS 
Hl FREQ: source frequency Hs >o 
COFF: integration contour offset - dB/A COFF > 0 
a a a A 7 PP SD SS SD 
IV NL: number of layers, incl. halfspaces - NL >2 
: depth of interface m - 
‘ cC: compressional speed m/s cc >0 
CS: shear speed m/s 
AC: compressional attenuation — dB/A ac >0 
“ AS: = shear attenuation dB/A as>0 
RO: density g/cm? no >0 
RG: rms value of interface roughness m - 
CL: correlation length of roughness m CL>0 
er EP 
Vv SD: source depth (mean for array) ™m - 
NS: number of sources in array - NS >0 
DS: vertical] source spacing m DS >0 
AN: grazing angle of beam deg - 
TA: = array type - is<macs 
: focal depth of beam m Fp # SD 
VI RDi: depth of first receiver m - 
RD2: depth of last receiver m RD2 > RD1 
wR: number of receivers” - NR >0O 
Ig: plot output increment - IR >0 
Vil CHIN: minimum phase velocity m/s CKIN > 0 
CHAX: maximum phase velocity m/s - 
Vint WW: number of wavenumber samples - NW = 2 
Ici: first sampling point - i>1 
+ Ic2: last sampling point - Ic2 < NW 





Figure 12 SAFARI Inputs [Ref. 7:p. 68] 


by its respective velocity and dividing by 1000), the density 
in g/cm’, and the RMS roughness which was left as zero for the 


study. For this study an average fluid solid interface depth 


of 3400 m was used in all runs. The next block contains the 





source depth which was 122 m for the study. The sixth block 
contains the receiver depth of the uppermost receiver, the 
lowermost receiver if more than one, the number of receivers 
and the interim receiver depths. Here all plots are of just 
one receiver depth and have the required default place holders 
of 1 for the number of receivers and interim receiver depths. 
The seventh block contains the minimum and maximum phase 
velocity speed. The eighth block contains the number of FFT 
sampling points, the starting sampling point, and ending 
sampling point. The ninth through thirteenth blocks deal with 
the various plotting parameters. Blocks nine, ten, and 
thirteen are used in generating the output data per the 
reference guide [Ref 7] and are shown in Figure 13. A 
representative file sample using contour offset is shown in 


Figure 14. 
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Block Parameter Units Limits 














IX RMIN: minimum range on plots km - 
RMAX: maximum range on plots km - 
RLEN: length of z-axis for all plots cm RLEN > 0 
RINC: distance between tick marks km RINC > 0 
x! TMIN: minimum transmission loss dB - 
TMAX: maximum transmission loss dB - 
TLEN: length of vertical TL axes cm TLEN > 0 
TINC: distance between tick marks dB TINC > 0 
xP DCUP: minimum depth for plots m - 
DCLO: maximum depth for plots m - 
DCLN: length of depth axis cm DCLN > 0 
DCIN: distance between tick marks m DCIN>0 
xIr ZXIN; minimum contour level dB - 
Z2MAX: maximum contour level dB - 
ZINC: contour level increment dB ZINC > 0 
xi1* VLEF: wave speed at left border m/s - 
VARIG: wave speed at right border m/s - 
VLEN: length of wave speed axis cm VLEN > 0 
VINC: wave speed tick mark distance m/s VINC > 0 
DVUP: depth at upper border m - 
DVLO: depth at lower border m - 
DVLN: length of depth axis cm DVLN > 0 
DVIN: depth-axis tick mark interval m DVIN > 0 





* Only for options 4, D, I and T. ? Only for options C and D. * Only for option C. 
* Only for option 2. 


Figure 13 SAFARI Plot Parameters (Ref. 7:p. 69] 
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Q 1544.88 ~1542.97 0 0 
2. ° 1542.97 -1543.93 0 


1543.93 -1542.00 0 
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-1540.93 0 
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1539.59 


1 


Figure 


538.72 


1530.84 
1530.14 
1530.56 
1528.08 
1527.04 
1526.09 
1526.05 
1523.25 
1519.53 
1515.17 
1511.40 
1503.42 


1496.45 -1494.06 


435.32 


4096 
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14 SAFARI Input File with standard 


-1537.45 0 
1537.45 ~1530.84 0 


-1530.14 
~1530.56 
~1528.08 
-1527.04 
-1526.09 
-1526.05 
-1523.25 
-1515.17 
~1511.40 
-1503.42 
~1496.45 


-1493.51 
-1497 .37 
-1497.44 
-1497.85 
-1500.13 
~1501.43 
~1506.97 
1511.75 
-1517.13 
~1524.38 
116.00 
122.00 
128.00 
134.00 
146.02 
158.02 
401.77 
552.79 
672.14 


oo0o0000 
ft te bet bt 
ooo0oo00o0o 


ooooqcooo0oo0ooo°ce 
qoooo0oooeoooc]e 
Fate bat ee at et ak pt ed et 
qooocooocoooco 


ceooocoo0oooecde 
oooooqocoocos 
Ped feed fed ad at bt bt pet fet bet bes 


ococoocoooeso 


0.2354 


1.5428 
1.6568 


0.0771 
0.0792 
0.0814 1.7766 
0.0837 1.9015 
0.0881 2.1567 
0.0926 2.4256 
0.1341 
0.2349 0 2.00 
0 2.18 


1.46 
1.47 
1.48 
1.49 
1.51 
1.53 


0 1.69 0 


0 
0 


1074.68 0.1353 0 2.32 0 
5175.0 2725.0 0.1035 0.19075 2.75 


11 


contour offset selected 
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V. MODEL INVESTIGATION 
In the SAFARI model, two issues must be investigated. The 
first is whether k,,, and k,,, the two endpoints of the sampling 
regime alias the result by truncating the depth-dependent 
Greens function at a point where the contribution from the 
integrand is significant. If the integration is truncated at 
a point where the kernel’s amplitude is significant the 
discontinuity will give rise to aliasing. The second is the 
whether a fine enough sampling of the horizontal wavenumber k 
is taking place, to ensure no aliasing is occurring in the 
depth-dependent Green’s function (the integrand), which 
oscillates sharply over k,,,-k,,, the span of the horizontal 
wavenumber. As previously stated, the contour offset is used 
to smooth the kernel by stepping off the axis. This has the 
effect of reducing the maximum spatial frequency at which 
aliasing commences, and therefore in range reduces the 
distance at which the onset of aliasing occurs. 
A. MODEL ALIASING 
1. Aliasing Background 
A standard definition of aliasing is 

"if a signal contains no frequency components above a 

frequency f,,, the signal can be uniquely represented by 

equally spaced samples if the sampling frequency f, is 


greater than twice f,,,. That is, the sampling frequency 
must satisfy the inequality f, > 2f,,,."[Ref. 12] 
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The inequality is known as the Nyquist rate. Since frequency 
is analogous to range and time analogous to wavenumber, 
frequency aliasing occurs when the time domain is 
undersampled. Similarly, in SAFARI, range aliasing occurs 
when wavenumber space is undersampled. Since the wavenumber 
sampling interval Ak=(k,,,-k,,,) /M, where M is the total number 
of sampling points, either increasing M or decreasing the 
wavenumber span k,,,-k,,, will reduce aliasing in range. Just 
as an impulse function in time domain yields a constant level 
in the frequency domain, so to will sharp peaks in wavenumber 
space yield a constant level in the spatial domain. If the 
Nyquist rate is violated, i.e., the frequency k, is not 
greater than twice k,,, the resulting transmission loss curve 
will show aliasing. 
2. Aliasing Results 

In SAFARI, the horizontal wavenumber k, governs the 
distance at which the aliasing commences. Given the rapidly 
varying integrand resulting from Eq. 3.27, the wavenumber 
range must be large enough and the sampling period fine enough 
to allow for accurate integration. There are five model input 
parameters governing the wavenumber integration: 

®c,, the minimum phase velocity 
® Cc... the maximum phase velocity 


e NW the number of wavenumber samples, which must be an 
integer power of 2, i.e. NW=2™, 


# 1ICl1 the first sampling point 





e1C2 the last sampling point 
The relationship between the limits of the horizontal 
wavenumbers k and phase speeds cC is Keg=2mf [Com and 
Knox=2UE£/Con, Where f is the frequency of interest. The 
wavenumber interval Ak is determined by taking (k,,,-k,,)/ M, 
where, M is the total number of sampling points. Phrased in 
another fashion, k,,, = Kum + (M-1)4k. Then recalling Eq 3.26 
Ar can be determined and the window range R =MAr can be 
determined. This is the maximum range at which the integrand 
can not be aliased due to the truncation. In this study, all 
graphs start with IC1 equal to 1 and IC2 equal to the last 
sampling point M. To ensure all compressional and shear 
speeds were incorporated, C,, was chosen well below the 
minimum speed in the profile. Unless there are computing 
constraints, a very large value should be chosen for C,,, to 
ensure truncation does not occur too soon. The first graphs 
Figure 15 and Figure 16 show the transmission loss curve and 
the wavenumber integrand. The transmission loss curve shows 
the beginning of aliasing just before 30 km, which supports 
the cutoff of the model at its predicted value of 29000 nm. 
This number was. arrived at via Equation 3.26. ,i.e., 
MAr=27/Ak, where Ak is k,,, - Knuae The wavenumber integrand 
shows a very confined area over which a very abeape series of 
oscillations take place on the axis of integration. Figure 


17 and Figure 18 compare the effect of reducing the number of 
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sampling points by half but keeping Ak the same by reducing 
Knax7Kmin iM half. There are two items to notice. First the 
onset of aliasing occurs at the predicted range of 29,000 n, 
the same as the previous eecuit: Second, the graphs are not 
the same, indicative of a reduction in aliasing in the kernel 
when the higher value of C,;,, was used, as is expected since 
the span in k has been reduced. 

Inspection of a series a series of these curves allows 
a C,,, and C,, which will provide the greatest alias-free range 
for a given number of sampling points. Figure 19 where 
Cnn=740 m/s, and C,,=1E8 m/s (1X10*), show no aliasing but the 
picture has converged to a more stable version. In this case, 
values of phase speed below 740 m/s and above 1E8 m/s do not 
significantly contribute. The integrand in Figure 20 has been 
rescaled automatically by SAFARI. There still is a large 
degree of ‘’peakness’ and this will be investigated later in 


the contour offset series. 


B. CONTOUR OFFSET 

| As previously described the contour offset is designed to 
smooth out the peaks in the integrand by moving away from the 
poles into the complex plane. This is analogous to going 
along the side of the mountain as opposed to going over its 
peaks. Since the integrand is composed of many spatial 


frequencies and the kernel is smoothed by taking an offset 


contour, this is the same as reducing the maximum range at 








which the onset of aliasing will occur. Given the number of 
sampling points is not changed, nor is the span of k, aliasing 
should begin at smaller values of r. Figure 21 and Figure 22 
show the onset of aliasing as moving to shorter ranges we the 
contour is moved off the axis. Figures 23 through 29 show the 
effect of using the contour offset for the 43 layer model with 
Crin=50 m/S and C,,=1E8 m/s. Figure 23 has no offset and each 
following Figure has an increase in the contour offset by a 
factor of 0.25 times the standard contour offset. With a 
standard offset and 4096 points, the aliasing begins at about 
25 km. Doubling the number of points to 8192 moves the 
aliasing out beyond the limits of the graph, as seen in Figure 
30 and Figure 31. 

In Figure 19 and Figure 20, C,,, was 1E8 and C,;,, was 740 
with no contour offset used. With a standard contour offset, 
Figure 32 shows the change in the nulls. This results ina 
small change in the transmission loss curve as shown in Figure 
33; the nulis are deeper and the overall shape shows a few 
more peaks. In the next series, the standard contour offset 
was used, and C,, was held constant at 740 m/s, but C,,, was 
reduced to 3400 m/s (Figure 34 and Figure 35). This had little 
effect on the results indicating that phase speeds above 3400 
m/s contribute little to the propagation. Figure 35 bears this 
out as small k shows no significant contribution to 


propagation. The largest effect was between 5 to 10 km. The 
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next comparison was in running the same profile with 1024 
points vice 4096. Figure 34 shows a transmission loss curve 
almost identical to the 4096 pt curve in Figure 34, with some 
very minor changes in the nulls around 8 km. The integrand in 
Figure 37 has changed dramatically. What were originally 
spikes in Figure 20, now have some width and continuity, 
making for an easier evaluation of the integrand and thus less 
aliasing due to the kernel. Next is the comparison of Figure 
36 with the 8192 point graph in Figure 31 both using the 
contour offset feature, to observe the losses due to the 
change in parameters. Clearly not much is lost in the overall 
shape of the transmission loss curve. A few more features are 
present with deeper nulls and smoother peaks and some 
resolution of a few peaks around 6 km, but otherwise the shape 
is the same, with the 1024 point graph taking less time to 
produce. Finally increasing the number of sampling points to 
16384 points (Figure 38) yielded no change in the graph over 
that of the 8192 point graph (Figure 31) over the range of 
interest. At this point the model has converged. Note the 
difference in the integrands between Figure 39 and Figure 37. 

Finally, Figure 40 shows the difference between the water 
column alone and using the water column and sediment combined. 
Clearly sediment contribution to the final transmission loss 


is important. 
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VI. COMPARISON OF SAFARI WITH SWALLOW FLOAT O DATA 

In this chapter the results of the Swallow Float records 
are presented graphically for 7 Hz and 10 Hz both without and 
with range offsets from regions 1 and 2. The results are 
visually compared and contrasted. Next the data is compared 
with the model results. In all cases Matlab 386 was used for 
the graphic presentations. All figures are located in 
Appendix D. 
a. 7 HZ DATA 

Figure 41 shows the 7 Hz transmission loss data taken 
during region 1. Figure 42 shows the 7 Hz transmission loss 
data taken during region 2. Figure 43 compares region 1 and 
region 2 data by offsetting region 1 to match predominant 
features. This offset was 2.4 km. 
B. 10 HZ DATA 

Figure 44 shows the 10 Hz transmission loss data taken 
from region 1. Figure 45 shows the 10 Hz transmission loss 
data taken from region 2.) Figure 46 compares region 1 with 
region 2 by offsetting region 1 to match predominant features. 
This offset was 2.5 km. 
C. SAFARI MODEL VS 7 HZ DATA 

Figure 47 shows the SAFARI model at 7 Hz vs 7 Hz 
transmission loss data taken from region 1. Figure 48 shows 


the SAFARI model at 7 Hz vs 7 Hz transmission loss data taken 
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from region 2. Figures 49 and 50 compare the SAFARI model at 
7 Hz vs 7 Hz transmission loss data taken from regions 1 and 
2 respectively, offset in range to match predominant features. 
The offset was 2.4 km for region 1 and none for region 2. 
D. SAFARI MODEL VS 10 HZ DATA 

Figure 51 shows the SAFARI model at 10 Hz vs 10 Hz 
transmission loss data taken from region 1. Figure 52 shows 
the SAFARI model at 10 Hz vs 10 Hz transmission loss data 
taken from region 2. Figures 53 and 54 compare the SAFARI 
model at 10 Hz vs 10 Hz transmission loss data taken from 
regions 1 and 2 respectively, offset in range to match 


predominant features. The offset was 0.5 km for region 1 and - 


0.75 km for region 2. 





VII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The initial findings in the Swallow float 0 data show 
excellent agreement with the SAFARI model; the model produces 
a result that compares well with ocean data. The model shows 
that there is a distinct problem with the range used for the 
data recorded in region 1. No explanation was found as to why 
the offset found in region 1 at 7 Hz was much larger than the 
navigation error mentioned in Reference 4. The data at 7 Hz 
in region 2 and at 10 Hz in both regions 1 and 2 match up 
reasonably well given errors in the navigation. The slope of 
the bottom appears to have very little effect on the 
transmission loss in this case. By comparing results obtained 
from modelling the velocity profile of the water column alone 
versus one with sediments included, the impact of the bottom 
and the importance of modelling the bottom become apparent. 

SAFARI is extremely sensitive to aliasing based on its 
input values of phase speed, number of sampling points, and 
limits of integration. The contour offset does provide a way 
to generate less aliasing from poles at the cost of aliasing 


‘at long ranges. This can be corrected by increasing the 


number of sampling points. 








B. RECOMMENDATIONS 
Two areas which can be addressed in future theses are 
investigation of data obtained on other floats and 


investigation of the effect the starting bottom depth has on 
the model. 
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APPENDIX A 


A. NAVIGATIONAL DATA 

The record of the position of the R/V Gyre is enclosed. 
It is uncorrected with respect to the navigational offset 
mentioned in Chapter 2. The last page has the time the 


changes in navigational source occurred. 
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APPENDIX B 


A. SOUND VELOCITY PROFILES 

Representative samples of the sound velocity profile of 
the water column used in the SAFARI model are presented in 
this appendix. A full listing of all the samples made is 


available in Reference 1. 
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SITE "D" PROVINCE PROFILES: 


SITE "D" PROVINCE ONE SITE "D" PROVINCE TWO 
31. 4833, mt wh 90 8 3850 40 em 31.5592 mie EEE * 8 3850 32 2 


1342, 99 
; ; ‘ 1541.80 
43.6 26.16 37.44 1540.91 
51.7 25.32 37.43 = 1539.12 






54.1 24.94 37.42 1538.28 
54.7 24.69 37.39 1537.68 
59.8 24.35 37.16 1536.68 
67.6 24.12 36.90 1535.97 
80.9 23.24 36.78 1533.94 
95.9 22.45 36.75 1532.17 
99.7 22.42 36.76 1532.16 
104.2 22.01 36.76 1531.21 
114.6 21.73 36.74 1530.64 
118.3. 21.40 36.74 1529.83 
134.3. 21.10 36.71 1529.28 
162.3 19.97 36.71 1526.71 

206.3 19.26 36.70 1525.50 

267.5 19.00 36.61 1525.66 





428.0 17.07 36.46 1522.55 
561.8 14.35 36.22 1516.06 
622.0 13.57 36.02 1514.26 
668.1 12.18 35.86 1510.16 
691.3 11.86 35.78 1509.35 
763.8 9.99 35.60 1503.68 
868.0 8.02 35.29 1497.69 





































900.0 7.75 35.23 1497.10 
913.0 7.41 35.19 1495.98 
989.0 6.50 35.10 1493.54 
1090.0 6.04 35.08 1493.37 
1160.0 5.64 35.07 1492.93 
1290.0 5.02 35.04 1492.54 1524.38 
1416.0 4.60 35.01 1492.88 
1628.0 4.19 34.99 1494.70 
2041.0 3.81 34.98 1500.13 
2176.0 3.58 34.96 1501.43 
2611.0 3.13 34.95 1506.97 
3018.0 2.62 34.92 1511.75 

0 2.33 1517.13 

2.18 1524.38 
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PROVINCE THREE PROFILE 
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PROVINCE FOUR PROFILE 
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APPENDIX C 


A. GEOACOUSTIC BOTTOM 

The geoacoustic bottom model and its results are presented 
in this appendix and are taken from Reference 1. The 
conversion of the compressional and shear attenuation 
coefficients from dB/(m*kHz) to dB/A, needed a an input in 
SAFARI, is as follows: 

A,=A,,*C,/1000 

where A, is the attenuation coefficient of the respective 
compressional wave or shear wave in the SAFARI model, A, is 
the attenuation coefficient of the geoacoustic model Cc, is the 


velocity of the x wave, where x is compressional or shear, and 


1000 is a scaling factor from kHz to Hz. 





Unconsolidated layer - Calcareous Sediment, (DSDP Well 533) 
Vp = 1.455 + 1.713() - 0.374(D)?_ (Hamilton, 1980) 


Vs = 110.0 + 1.2005(@D) (for Vp < 1.512) (Hamilton, 1976) 

Vs = (3.884 x Vp) - 5.757 (for Vp 1.512-1.555 km/s) (Hamilton, 1980) 
Vs = 0.991 - (1.136 x Vp) + 0.47(Vp)* (for Vp 1.65-2.15 km/s) 

Vs = (0.78 x Vp) - 0.962 (for Vp > 2.15 km/s) 


Kp: Estimated from Hamilton, 1980 
Ks: Estimated from Hamilton, 1980 


p = 1.1379 (Vp) - 0.20484 (Hamilton, 1978) 


Gas layer beneath BSR (m/s)- 
Vp = 1090.4 + 1.2892(D) - [(2.8628 x10 )x(D)?] 


Vs = Vp/1.9 (Hamilton, 1980; Pandit and King, 1982) 


Kp: Estimated from Hamilton, 1980 
Ks: Estimated from Hamilton, 1980 


p =0.112 (Vp) + 2.088 (Hamilton, 1978) 


BWSS = Bottom Water Sound Speed 
D=Depth 
Vp = Compressional Velocity 
Vs = Shear Velocity 
= Density 
z Compressional Attenuation 
Ks = Shear Attenuation 
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Geoacoustic models for Site “D" 
Event 1Leg 1 Start 31°26'13" N, 74°49'27" W 
Bearing 193° Range 0 nm 
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Ks: 0.07 dB/m*kHz ; 
p: 2.75 g/cc 
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Appendix G. Event ILeg 1Start 31°26 '13" N, 74°49 '27" W 
‘Bearing 193° ‘ Range 10.45 nm 



















(m) (m/s) (m/s) (dB/m*kHz) (dB/m*kHz) g/cc 
Q 455.00 110.00 0.0 OO 1.4 
0 1463.56 116.00 0.0527 13.30 1.46 

10.0 1472.09 122.00 0.0538 13.58 1.47 

15.0 1480.61 128.00 0.0550 13.88 1.48 

20.0 1489.11 134.00 0.0562 14.19 1.49 

30.0 1506.05 146.02. 0.0585 14.77 1.51 

40.0 1522.92 158.02 0.0608 15.35 1.53 

50.0 1539.72 223.27 0.0632 15.95 1.55 

75.0 1581.37 313.02 0.0690 17.42 1.59 

100.0 1622.56 359.85 0.0748 18.88 1.64 
125.0 1663.28 401.77 0.0806 20.00 1.69 
150.0 1703.54 419.74 0.0865 20.00 1.73 
200.0 1782.64 459.49 0.0981 20.00 1.82 
250.0 1859.88 503.98 0.1098 20.00 1.91 
300.0 1935.24 552.79 0.1214 20.00 2.00 
400.0 2080.36 661.82 0.1575 20.00 2.10 
500.0 2218.00 768.04 0.1425 20.00 2.17 
600.0 1761.12 434.78 0.2550 30.00 1.66 
625.0 1910.65 440.69 0.1875 25.00 1.90 _ 
650.0 2060.11 645.42 0.1200 20.00 2.14 
700.0 2092.86 672.14 0.1125 20.00 2.18 
800.0 2155.54 719.32 0.0975 20.00 2.25 
900.0 2215.74 766.28 0.0825 20.00 2.26 
1000.0 2275.00 812.50 17.04 2.27 
1500.0 2611.13 1074.68 0.05181 3.08 2.32 
3155.00 1498.90 - 9.11 2.37 
3836.13 2030.18 5.55 2.42 





Basement: 
Vp: 5175 m/s 
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Appendix G. Event Leg 1Start 31926 '13" N, 74°49 '27" W 
Bearing 193° Range 14.21 nm 


(m) (m/s) (m/s) (dB/m*kHz) (dB/m*kHz) g 
( 5.0 0.0 6.051 13.00 


110.00 


g°> 












. 0 "4 
0 1463.56 116.00 0.0527 13.30 1.46 
10.0 1472.09 122.00 0.0538 13.58 1.47 
15.0 1480.61 128.00 0.0550 13.88 1.48 
20.0 1489.11 134.00 0.0562 14.19 1.49 
30.0 1506.05 146.62. 0.0585 14.77 1.51 
40.0 1522.92 158.02 0.0608 15.35 1.53 
50.0 1539.72 223.27 0.0632 15.95 1.55 
75.0 1581.37 313.02 0.0690 17.42 1.59 ; 
100.0 1622.56 359.85 0.0748 18.88 1.64 ‘ 
125.0 1663.28 401.77 0.0806 20.00 1.69 
150.0 1703.54 419.74 0.0865 20.00 1.73 
200.0 1782.64 459.49 0.0981 20.00 1.82 » 
250.0 1859.88 503.98 0.1098 20.00 1.91 
300.0 1935.24 $52.79 0.1214 20.00 2.00 
400.0 1560.27 330.91 0.3150 30.00 1.58 
425.0 1659.06 389.67 0.2120 25.00 1.69 
500.0 1953.88 565.68 0.1425 20.00 2.02 
580.0 2012.19 608.14 0.1305 20.00 2.08 
600.0 2026.15 618.78 0.1275 20.00 2.10 
650.0 2060.11 645.42 0.1200 20.00 2.14 
700.0 2092.86 672.14 0.1125 20.00 2.18 
800.0 2155.54 719.32 0.0975 20.00 2.25 
900.0 2215.74 166.28 0.0825 20.00 2.26 
1000.0 2275.00 812.50 0.0675 17.04 2.27 
1500.0 2611.13 1074.68 0.0518 13.08 2.32 
2000.0 3155.00 1498.90 0.0361 . 2.37 
3836.13 . 2.42 









2030.18 0.0220 





Appendix G. Event 1Leg 1Start 31°26 '13" N, 74°49 27" W 









“Bearing 193° Range 28.63 nm 
(m) (m/s) (m/s) (dB/m*kHz) _(dB/m*kHz) g/cc 
5.0 455.00 0.00 0.051 13.00 1.45 
5.0 1463.56 116.00 0.0527 13.30 1.46 
10.0 1472.09 122.00 0.0538 13.58 1.47 
15.0 1480.61 128.00 0.0550 13.88 1.48 
20.0 1489.11 134.00 0.0562 14,19 1.49 
30.0 1506.05 146.02. 0.0585 14.77 1.51 
40.0 1522.92 158.02 0.0608 15.35 1.53 
50.0 1539.72 223.27 0.0632 15.95 1.55 
75.0 1581.37 313.02 0.0690 17.42 1.59 
2 100.0 1622.56 359.85 0.0748 18.88 1.64 
125.0 1663.28 401.77 0.0806 20.00 1.69 
. 150.0 1703.54 419.74 0.0865 20.00 1.73 
200.0 1782.64 459.49 0.0981 20.00 1.82 
250.0 1859.88 503.98 0.1098 20.00 1.91 
300.0 1935.24 $52.79 0.1214 20.00 2.00 
400.0 2080.36 661.82 0.1575 20.00 2.10 
500.0 2218.00 768.04 0.1425 20.00 2.17 
600.0 2348.16 869.56 0.1275 20.00 2.21 | 
700.0 2470.84 965.26 0.1125 20.00 2.24 
800.0 2586.04 1055.11 0.0975 20.00 2.26 
900.0 2693.76 1139.13 0.0825 20.00 2.28 
1000.0 2794.00 1217.32 0.0675 17.04 2.29 
1500.0 3183.00 1520.74 0.0518 13.08 2.34 
2000.0 3385.00 1678.30 0.0361 9.11 2.36 
2500.0 3955.00 2122.90 0.0204 5.15 2.42 
2930.0 4385.00 2458.30 0.0070 2.47 








Basement: 
p: 5175 m/s 
Vs: 2725 m/s 


Kp: 0.02 dB/m*kHz © 
Ks: 0.07 dB/m*kHz 
p: 2.75 g/cc 





Appendix G. Event 1Leg 1Start 31926 '13" N, 74°49 '27" W 


Bearing 193°Range 40.13 nm 
D Vp Vs p Ks p 
(m) (m/s) (m/s) (dB/m*kHz) _(dB/m*kHz) g/cc 
0.0 455.00 10.00 0.051 13.00 1.4 
5.0 1463.56 116.00 0.0527 13.30 1.46 
10.0 1472.09 122.00 0.0538 13.58 1.47 
15.0 1480.61 128.00 0.0550 13.88 1.48 
20.0 1489.11 134.00 0.0562 14.19 1.49 
30.0 1506.05 146.02. 0.0585 14.77 1.51 
40.0 1522.92 158.02 0.0608 15.35 1.53 
50.0 1539.72 223.27 0.0632 15.95 1.55 
75.0 1581.37 313.02 0.0690 17.42 1.59 
100.0 1622.56 359.85 0.0748 18.88 1.64 e 
125.0 1663.28 401.77 0.0806 20.00 1.69 j 
150.0 1703.54 419.74 0.0865 20.00 1.73 # 
2 1782.64 459.49 0.0981 20.00 1.82 6 
250.0 1859.88 503.98 0.1098 20.00 1.91 ie 
300.0 1935.24 552.79 0.1214 20.00 2.00 f 
400.0 2080.36 661.82 0.1575 20.00 2.10 . 
500.0 2218.00 768.04 0.1425 20.00 2.17 x 
580.0 1741.57 424.64 0.2610 30.00 1.65 ‘2 
605.0 1952.21 503.44 0.2110 25.00 1.82 & 
650.0 2060.11 645.42 0.1200 20.00 2.14 ae 
700.0 2092.86 672.14 0.1125 20.00 2.18 Si 
800.0 2155.54 719,32 0.0975 20.00 2.25 = 
900.0 2215.74 766.28 0.0825 20.00 2.26 i 
1000.0 2275.00 812.50 0.0675 17.04 2.27 rf 
1500.0 2611.13 1074.68 0.0518 13.08 2.32 i 
2000.0 3155.00 1498.90 0.0361 9.11 2.37 : 
2450.0 3836.13 2030.18 0.0220 5.55 2.42 ‘ 
2500.0 3861.13 2049.68 0.0205 5.17 2.43 ie 
2675.0 3948.63 2117.93 0.0150 3.79 2.44 i 
2700.0 3961.13 2127.68 0.0142 3.58 2.45 y 
2800.0 4011.13 2166.68 0.0111 2.80 2.46 3 
3000.0 4111.13 2244.68 0.0048 1.21 2.48 us 
Basement: es 
Vp: 5175 mvs ie 
Vs: 2725 ms é 
Kp: 0.02 dB/m*kHz 


Ks: 0.07 dB/m*kHz . eo 
9: 2.75 g/cc oe 





Appendix G. Event 1Leg 1Start 31°26 '13" N, 74°49 '27" W 
- ‘Bearing 193° Range 42 nm 
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APPENDIX D 


A. MODEL RESULTS 

Figures 15 through 40 show the results of the SAFARI model 
run at a frequency of 7 Hz. OASTL STLRAN and OASTL SINTGR 
stand for Ocean Acoustic and Seismic Exploration Synthesis 
Transmission Loss, Single frequency Transmission Loss versus 
Range, and Ocean Acoustic and Seismic Exploration Synthesis 
Transmission Loss, Single frequency Integrand. They are 


located on the top right part of the graphs. 
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Figure 15 Transmission Loss SAFARI Float 0 
Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 mn, 
Max FFT range=29.250 km 
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Figure 16 Integrand SAFARI Float 0 
Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 nm, 
Max FFT range= 29.250 km 
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Figure 17 Transmission Loss SAFARI Float 0 
Cmin=100 m/s, Cmax=2E8 m/s, 2048 points, Ar=14.279 m, 
Max FFT range=29.243 km 
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Figure 18 Integrand SAFARI Float 0 
Cmin=100 m/s, Cmax=2E8, 2048 points, Ar=14.279 m, 
Max FFT range=29.243 km 
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Figure 19 Transmission Loss SAFARI Float 0 
Cmin=740 m/s, Cmax=1E8 m/s, 4096 points, Ar=105.7 nm, 
Max FFT range=432.903 km 
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Figure 20 Integrand SAFARI Float 0 
Cmin=740 m/s, Cmax=1E8 m/s, 4096 points, Ar=105.7 nm, 
Max FFT range=432.903 km 
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Figure 21 Transmission Loss SAFARI Float 0 

No Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 m, 
Max FFT range=29.250 km 


77 





OASTL STLRAN 


Tp) 
Y) 
3 
Za 
O 
Y) 
WY 
a 
YD 
Zz 
Ll 
ag 
|— 


SAFARI FLOAT 0 


PdL/AP) SSe.}S JOLUION 





Figure 22 Transmission Loss SAFARI Float 0 

0.25 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 nm, 
Max FFT range=29.250 km 
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Figure 23 Transmission Loss SAFARI Float 0 

No Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 mn, 
Max FFT range=29.250 km 
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Figure 24 Transmission Loss SAFARI Float 0 

0.25 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 m, 
Max FFT range=29.250 km 
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Figure 25 Transmission Loss SAFARI Float 0 

0.5 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 m, 
Max FFT range=29.250 km 
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Figure 26 Transmission Loss SAFARI Float 0 

0.75 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 m, 
Max FFT range=29.250 km 
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Figure 27 Transmission Loss SAFARI Float 0 
1.0 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 m, 
Max FFT range=29.250 km 
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Figure 28 Transmission Loss SAFARI Float 0 
1.25 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 n, 
Max FFT range=29.250 km 
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Figure 29 Transmission Loss SAFARI Float 0 
1.50 Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 mn, 
Max FFT range=29.250 km 
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Figure 30 Transmission Loss SAFARI Float 0 
Standard Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 n, 
Max FFT range=29.250. km 


86 








OASTL STLRAN 


Cn) 
Y) 
G 
eae 
O 
YY) 
Y) 
> 
D 
Ps 
— 
i 
}— 


SAFARI FLOAT O 


(Odi //GP) SS24}S JOLLUON 





Figure 31 Transmission Loss SAFARI Float 0 


Standard Contour Offset 
Cmin=50 m/s, Cmax=1E8 m/s, 8192 points, Ar=7.141 m, 


Max FFT range=58.500 km 
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Figure 32 Integrand SAFARI Float 0 

Standard Contour Offset 

Cmin=740 m/s, Cmax=1E8 m/s, 4096 points, Ar=105.7 mn, 
Max FFT range=432.903 km 
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Figure 33 Transmission Loss SAFARI Float 0 
Standard Contour Offset 

Cmin=740 m/s, Cmax=1E8 m/s, 4096 points, Ar=105.7 m, 
Max FFT range=432.903 km 
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Figure 34 Transmission Loss SAFARI Float 0 

Standard Contour Offset 

Cmin=740 m/s, Cmax=3400 m/s, 4096 points, Ar=105.7 m, 
Max FFT range=432.903 km 
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Figure 35 Integrand SAFARI Float 0 

Standard Contour Offset 

Cmin=740 m/s, Cmax=3400 m/s, 4096 points, Ar=105.7 m, 
Max FFT range=432.903 km 
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Figure 36 Transmission Loss SAFARI Float 0 

Standard Contour Offset 

Cmin=740 m/s, Cmax=3400 m/s, 1024 points, Ar=105.7 m, 
Max FFT range=432.903 km 
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Figure 37 Integrand SAFARI Float 0 

Standard Contour Offset 

Cmin=740 m/s, Cmax=3400 m/s, 1024 points, Ar=105.7 mn, 
Max FFT range=432.903 km 
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Figure 38 Transmission Loss SAFARI Float 0 

Standard Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 16384 points, Ar=7.141 nm, 
Max FFT range=117.000 km 
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Figure 39 Integrand SAFARI Float 0 

Standard Contour Offset 

Cmin=50 m/s, Cmax=1E8 m/s, 16384 points, Ar=7.141 n, 
Max FFT range=117.000 km 
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Figure 40 Transmission Loss SAFARI Float 0 

Water Column Only 

Cmin=50 m/s, Cmax=1E8 m/s, 4096 points, Ar=7.141 m, 
Max FFT range=29.250 km 
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Figure 41 Transmission Loss, NATIVE 1, Event 1, 
7 Hz, Region 1 
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Figure 42 Transmission Loss, NATIVE 1, Event 1, Float 0, 
7 Hz, Region 2 
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Figure 43 Transmission Loss, NATIVE 1, Event 1, Float 0, 
7 Hz, Region 1 vs Region 2 
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Figure 44 Transmission Loss, NATIVE 1, Event 1, Float 0, 
10 Hz, Region 1 
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Figure 45 Transmission Loss, NATIVE 1, Event 1, Float 0, 
10 Hz, Region 2 





Range (Km) 


Region 1-dashed line 


- 
2 
bs 
fo) 
iv 
5) 
2 
5) 
= 
{o} 
cI 
a 
be 
a 


Region 2-solid line 





SSO7] UOISSTUISUBL 





Figure 46 Transmission Loss, NATIVE 1, Event 1, Float 0, 


10 Hz, Region 1 vs Region 2 
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Figure 47 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 1 at 7 Hz 
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Figure 48 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 2 at 7 Hz 
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Figure 49 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 1 at 7 Hz, with 2.4 km Offset in Range 
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Figure 50 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 2 at 7 Hz, with No Offset in Range 
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Figure 51 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 1 at 10 Hz 
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Figure 52 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 2 at 10 Hz 
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Figure 53 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 1 at 10 Hz, with 0.5 km Offset in Range 
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Figure 54 Transmission Loss, SAFARI vs NATIVE 1, Event 1, 
Float 0, Region 2 at 10 Hz, with -0.75 km Offset in Range 
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